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^ : ABSTRACT 

I In principle, the tensor metric (gravity-wave) perturbations that arise in in- 

^ . flationary models can, beyond probing the underlying inflationary model, 

I provide information about the Universe: ionization history, presence of a 

cosmological constant, and epoch of matter-radiation equality. Because ten- 
sor perturbations give rise to anisotropy of the cosmic background radiation 
(CBR) solely through the Sachs- Wolfe effect we are able to calculate an- 
alytically the variance of the multipole moments of this part of the CBR 
temperature anisotropy. In so doing, we carefully take account of the contri- 
bution of tensor perturbations that entered the Hubble radius during both 
the matter-dominated and radiation-dominated epoch by means of a transfer 



function. (Previously, only those modes that entered during the matter era 
were properly taken into account.) The striking feature in the spectrum of 
multipole amplitudes is a dramatic fall off for I > vT+^lss, where 2;lss is 
the red shift of the last-scattering surface, which depends upon the ioniza- 
tion history of the Universe. Finally, using our transfer function we provide a 
more precise formula for the energy density in stochastic gravitational waves 
from inflation, and, using the Cosmic Background Explorer Differential Mi- 
crowave Radiometer (COBE DMR) quadrupole normalization, we express 
the this energy density in terms of the "tilt" of the spectrum of tensor per- 
turbations alone and show that it is unlikely that the stochastic background 
of gravity waves can be detected directly in the foreseeable future. 



1 Introduction 



Quantum fluctuations arising during inflation lead to a spectrum of scalar 
(density) 0] and tensor (gravity-wave) [0] metric perturbations which are 
nearly scale invariant 0. In turn, both give rise to temperature anisotropy 
in the cosmic background radiation (CBR), with their relative contributions 
depending upon the steepness of the inflationary potential |Q, |^. CBR 
anisotropy and other astrophysical data, e.g., red-shift surveys, peculiar- 
velocity measurements, and data from gravity-wave detectors, can, in prin- 
ciple, be used to learn much about the inflationary potential in the narrow 
interval that governs the modes that affect astrophysically interesting scales. 
For example, they can be used to infer the value of the inflationary potential, 
its steepness and the change in its steepness 0, ^ . 

In all likelihood CBR anisotropy provides the cleanest and most sensitive 
probe of metric perturbations. The CBR anisotropics that arise due to scalar 
and gravity-wave fluctuations add incoherently and can thus be computed 
independently. The calculation of the anisotropy that arises due to scalar 
perturbations is both complicated and well understood: anisotropy arises 
due to at least three physical effects, the Sachs- Wolfe effect [0 (gravitational 
potential differences on the last-scattering surface), the velocity of the last- 
scattering surface, and intrinsic fluctuations in the CBR temperature at last 
scattering. Further, the ionization history and baryon density are very im- 
portant p|. (On large-angular scales, for standard recombination 9 > 2°, 
the Sachs- Wolfe effect dominates, and on small-angular scales the other two 
effects dominate.) 

The CBR anisotropy due to tensor perturbations arises solely from the 
Sachs- Wolfe effect. It depends significantly upon the red shift of the last- 
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scattering surface, and less significantly upon a possible cosmological con- 
stant and the red shift of matter-radiation equality (through the value of the 
Hubble constant). Thus, if the tensor contribution to CBR anisotropy can 
be separated from the scalar contribution, it provides a very direct probe 
of cosmology. The tensor contribution has been computed analytically, but 
only on large-angular scales, with other simplifying assumptions being made 
[|], and recently, it has been computed numerically in the case of standard 
recombination and a matter-dominated Universe . 



The purpose of our paper to give simple and accurate analytic formulae 
for the tensor contribution to the variance in the CBR temperature multi- 
poles (la^mP) ("angular power spectrum"). In the previous analytical work 
[^], only the modes that enter the horizon after matter-radiation equality 
were properly taken into account, so that these results are accurate only on 
large angular scales; we take into account the modes that enter the hori- 
zon before matter-radiation equality by means of a transfer function and 
thereby accurately describe the CBR anisotropy that arises on all angular 
scales. The transfer function also allows us to give a more precise expres- 
sion for the energy density in stochastic gravity waves produced by inflation, 
and, unfortunately, we show that it is very unlikely that this background 
can be detected directly. We compare our results for CBR anisotropy, where 
possible, to the numerical results of Ref. 0], and discuss how the tensor 
multipoles depend upon the underlying cosmological parameters: the red 
shift of last scattering, -Zlss! the red shift of matter-radiation equality, z^q] 
and the presence of a cosmological constant. 
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2 Gravity Waves and CBR Anisotropy 
2.1 Qualitative view 

In this Section we begin simply and gradually add detail, ending with our 
most general formulae. In that spirit we will assume scale-invariant metric 
perturbations to begin, and then generalize our results to allow for devia- 
tion from scale invariancc. For scale-invariant perturbations the amplitude 
of density perturbations and gravity- wave perturbations are independent of 
scale at horizon crossing in the post-inflationary Universe: ((5p/p)hor = £5 
and /iGW — £t- Horizon crossing is defined by /Cphys — k/ R H] R is the 
cosmic-scale factor and H is the expansion rate. Throughout we take the 
scale factor to be unity at the present epoch, so that comoving wavenumber 
k and physical wavenumber /cphys are equal today; comoving wavelength A 
and wavenumber are related by: \ = 2'K/k. 

It is useful to define the conformal time today, tq, and at the last- 
scattering event, tlss^ 



where we assume a flat, matter-dominated Universe today. The quantities 
To ~ 6000/i~^ Mpc and tlss — 6000/i~-^ Mpc/Vl + ^lss correspond to the co- 
moving size of the present horizon and that at last scattering. The comoving 
distance to the last-scattering surface c^lss = To(1 — 1/vT+^lss) — Tq, and 
thus the angle subtended by the horizon scale at last scattering corresponds 
to 6'lss ~ t-lss/tq (~ 2° for 2;lss = 1100). 

In the standard picture recombination occurs at a red shift of order 1300, 




(1) 
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i.e., the ionization fraction becomes small, and last-scattering occurs a red 
shift of about zlss — 1100, i.e., the photon mean-free path becomes greater 



than the Hubble scale |[Tl[. If the Universe remains ionized much later, 
or is reionized after recombination, the last-scattering event can occur much 
later, I + ^lss — (O.OSXefi^^)"^''^, where fi^ is the fraction of critical density 
contributed by baryons, Hq = lOO/ikms^^ sec^^ Mpc^^, and VIq = 1. Taking 
h = 0.4, Xe = 1, and saturating the primordial nucleosynthesis bound to the 
baryon mass density, fi^/i^ < 0.02 |12|, last scattering could occur as late 
as zlss — 76. While the conventional wisdom has it that cold dark matter 
models lack a plausible mechanism for reionizing, or keeping the Universe 
ionized at z < 1300, it has been suggested that a very early generation of 



stars could have reionized the Universe at red shift of order 100 or so |T3|. 

The physics underlying the Sachs- Wolfe effect is very simple: The tem- 
perature fluctuation on a given angular scale is roughly equal to the metric 
fluctuation on the corresponding length scale on the last-scattering surface. 
For tensor perturbations, the metric perturbation is equal to the gravity-wave 
amplitude. For scalar perturbations the metric perturbation is given by the 
fluctuation in the Newtonian potential; on length scales that have yet to 
enter the horizon at last scattering, or that entered the horizon after matter- 
radiation equality, the fluctuation in the Newtonian potential is given by the 
horizon-crossing amplitude of the density perturbation. For scales that enter 
the horizon while the Universe is still radiation dominated the fluctuation in 
the Newtonian potential decreases after horizon crossing, as R~^. 

The CBR temperature fluctuation that arises on a given angular scale 
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due to scalar perturbations through the Sachs- Wolfe effect only is roughly: 

\es e< Oeq; (3) 



'5T\ f5p\ 



^ J \ P J HOR,k{e) 

ST] ^ (k^) (SjA ^ M ^ 



where k{9) corresponds to the wavelength that subtends an angle 6 on the 
last-scattering surface, 

m - m .„-'. (4) 





Here, /ceq is the scale that crossed the horizon at matter-radiation equality, 
and 9eq ~ 1/vTT^q ~ 0.3°. 

For scale- invariant perturbations and 6 > 6eq ~ 0.3°, the Sachs- Wolfe 
contribution to the temperature fluctuation is independent of angular scale. 
On small-angular scales, 6 < ^eq, the Sachs- Wolfe contribution to the CBR 
temperature anisotropy decreases, but is a subdominant contribution to total 
anisotropy produced by scalar perturbations. On angular scales larger than 



the horizon at last scattering, 9 > ^lss ~ l/Vl + -^^lss (~ 2° for 2;lss = 1100), 
the Sachs- Wolfe effect is the dominant contribution to CBR anisotropy. 

The CBR temperature fluctuation on a given angular scale due to tensor 
perturbations arises through the Sachs- Wolfe effect and is roughly: 

(^^j ~/iGwkss,A;(^)]; (5) 

where hGw[zLSS,k{9)] is the (dimensionless) gravity- wave amplitude on the 
scale k{9) at last scattering. For gravity-wave modes that have yet to re- 
enter the horizon at last scattering, /igw(-2^lss, k) is just equal to ex', once a 
mode enters the horizon its amphtude red shifts as the scale factor, so that 
/iGw(-2^LSS) fc) = ^t(-Rhor/-Rlss)- For modes that enter the horizon during 



the matter-dominated epoch, -Rhor/-Rlss — (^tlss)~^; for modes that enter 
during the radiation-dominated epoch, -Rhor/-Rlss — {k / kEq) {kri^ss)'"^ ■ [A 
simple way to understand why the amphtude of a gravity-wave mode red 
shifts as the scale factor after it enters the horizon is that gravity-wave per- 
turbations correspond to massless bosons (gravitons) whose energy density 
(oc 'Tipi^fcpj^yg/iQ-^) red shifts as R""^.] 

For scale-invariant tensor perturbations the CBR anisotropy due to grav- 
ity waves is only independent of scale for 9 > O^ss', on smaller angular scales 
it decreases. 



On angular scales greater than that subtended by the horizon at last 
scattering, ^lss ~ 2° in the standard scenario, the ratio of the tensor to the 
scalar contributions to the temperature anisotropy is constant and equal to 
Et/ss] in turn, this ratio is related to the steepness of the inflationary poten- 
tial, evaluated when these scales crossed outside the horizon during inflation 
(about 50 e-folds before the end of inflation): St/es ~ 2:50 = {jn-piV /V)^q 
[^]. On smaller angular scales the tensor contribution falls as 9"^ because these 
scales are dominated by gravity-wave modes that have entered the horizon 
before the epoch of last scattering and have had their amplitudes red shifted. 
On the smallest angular scales the tensor contribution only decreases as 9, 
as these scales are dominated by gravity-wave modes that enter the horizon 
before matter radiation equality. The steep fall off of gravity-wave modes for 




(6) 
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I ^ V 1 + ^Lss provides a potential signature of 2;lss and tensor perturbations. 

It is conventional to expand the CBR temperature fluctuation on the sky 
in spherical harmonics: 

^ = E E aUr)Yimi^y, (7) 

where the unobservable monopole term and the dipole term, which is dom- 
inated by the contribution of the observer's peculiar velocity, are omitted. 
The multipole amplitudes depend upon the observer's position r. The quan- 
tity la^mp averaged over all observation positions (the ensemble average) is 
referred to as the angular power spectrum,]] and is related approximately to 
the CBR temperature fluctuation by 

2 

I ~/'(|a,„|2) for/- 20070. (8) 



ST' 



To be more precise, the rms temperature fluctuation averaged over the 
sky for a given experiment is given by 

where Fi is the appropriate response function for the experiment. For a two- 
beam experiment, where the temperature difference between two antennas 
of gaussian beam width a separated by angle 6 is measured. 



Fi = 2[l-Pi{cosd)]e 



-(«+l/2)V2 



For scale-invariant density perturbations the Sachs- Wolfe contribution to 
P times the angular power spectrum is roughly constant and equal to Ss 
^More precisely, the average is over all realizations of the fluctuation spectrum; we have 
implicitly assumed spatial ergodicity. 
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for / < Ieq] thereafter it decreases as For scale-invariant gravity-wave 
perturbations times the angular power spectrum is constant for I < /lss 
(~ 30 for standard recombination); it decreases as for / < Zeq; and and 
as Z"^ for / > /eq- 

Finally, if the scalar and tensor perturbations are "tilted" — that is not 
scale invariant — say es oc A"'^ and Et oc A"^, the previous results for {6T/T)g 
are modified by factors of 9°"^ and 6*°^ respectively, and for (|a;mP), by fac- 
tors of /"^""s g^Yid /~2ot respectively. The quantities as and ax are related 
to the power-law indices often used to characterize the scalar and tensor 
perturbations (see below) by 



The qualitative behaviour of the CBR anisotropy due to scalar and tensor 
perturbations is shown in Fig. 1. 

2.2 Quantitative view 

Some of what follows is a quick review of previous treatments included for 
completeness; for more details see Refs. P|. To begin, we write the line ele- 
ment for a flat, Friedmann- Robertson- Walker (FRW) cosmology in conformal 
form plus a small perturbation hfj_^: 



where rj^i, = diag(l, —1,-1,-1) and r is conformal time. Here we are only 
interested in gravity-wave perturbations and work in the transverse-traceless 
gauge, where the two independent polarization states are x and -|-, and 



as = {l-n)/2- 



ar = —nT/2. 



(10) 



It is simple to solve for the evolution of the cosmic-scale factor in terms 
of conformal time: 

R{t) = [t/to + i^Eq] ' - ^eq; (11) 

where we have assumed a flat Universe with matter and radiation (i.e., pho- 
tons with present temperature 2.73K and three masslcss neutrino species), 



i?o = 1, To = 2Hq^^1 + Req ~ 2Ho \ Req = 4.18 x lO^^ h''^ is the value 
of the cosmic scale factor at the epoch of matter-radiation equality, and 
teq — [V^ — Iji^EQ^o- (The conformal age of the Universe today differs from 
2Hq^ by a small amount due to the tiny contribution of the radiation energy 
density today: 



Ttoday — To 



+ - ^/R^ ^ 2H^' (l - yfR^) « 2H^' [1 - 0{1%)] , 



which we shall henceforth neglect.) 

We expand the gravity-wave perturbation in plane waves 

/i,,(x, r) = (27r)-^ J d'k h\,{T)e], e"^*^'- (12) 

where e*^. is the polarization tensor and i — x, -|-. The gravity- wave pertur- 
bation satisfies the massless Klein-Gordon equation 

K + 2 (I) K + k-'hi = 0; (13) 

where an overdot indicates a derivative with respect to conformal time and 
= k • k. 

The growing-mode solutions to this equation have simple qualitative be- 
haviour: before horizon crossing [kr <C 1) h\J{r) is constant; well after hori- 
zon crossing {kr ^ 1) oc cos kr/R. For modes that cross inside the 
horizon during the radiation-dominated era the exact solution prior to and 
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including the radiation-dominated era is jo{kT); for modes that cross inside 
the horizon during the matter-dominated era the exact solution is 3ji (kr) / kr. 
Here jo{z) = sin z/z and ji{z) = smz/z^ — cosz/z are the spherical-Bessel 
functions of order zero and one respectively, and both Bessel-function solu- 
tions have been normalized to unity for r ^ 0. 

Well into the matter-dominated era the temporal behaviour of modes 
that entered the horizon during the radiation-dominated era is also given by 
3ji{kT)/kT. Thus, for r ^ Teq the temporal behaviour of all modes is given 
by 3ji{kT)/kT, and it is useful to write 

hUr) = hm T{k/k^ci) (^^) ; (14) 

where the "transfer function" for gravitational waves, T{k/k^Q), is only a 
function of /c/Zceq, and 

kEQ = = (^1^ ^ 6.22 X lO-'h' Mpc-^ (15) 

is the scale that entered the horizon at matter-radiation equality. Note, 
during the oscillatory phase {kr ^ 1), 3ji(A;r)/A;r 3 cos kr / {krY ; in 
defining and computing the transfer function we have neglected the phase of 
the graviton oscillations, which, for our purposes, is not important. 



The transfer function has been calculated by integrating Eq. ([T3|) numer- 
ically from r = to r = To; a good fit to the transfer function is 

T(y) = [1.0 + 1.3% + 2.50?/2]V^ (16) 

where y = k/k-EQ- It might seem that one could have computed the trans- 
fer function at any time after the Universe becomes matter dominated and 
obtained the same result. However, as we shall emphasize again later, the 
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Universe becomes matter dominated more slowly than one might have ex- 
pected, and for this reason the transfer function calculated at an earlier epoch 
is different. Only well into the matter-dominated epoch, when the radiation 
content is very neghgible, does it take this functional form. Using the fact 
that once a mode is well inside the horizon /^^(t) oc cos kr/R it follows that 
for modes with k ^ /ceq the transfer function at an earlier epoch is related 
to that today by 

Tr{k/kEQ) = {r/rofR-'T{k/kEQ) = [l + 2x - 2^^+^] %{k/k^ci)] (17) 

where x = Req/R{t), is the transfer function at conformal time r, and 
T is the transfer function today. The evolution of the transfer function is 
shown in Fig. 2. 

Once a mode has crossed inside the horizon one can sensibly talk about 
the energy density in gravitons associated with that mode; it is given by 

21,5 



dpi _ mpi k 
"dk " 327r3i?2 l^^^^^l ' ^^^^ 



where \h\.\^ is the average of \h\.\^ over all directions k and over several 
periods. 

Inflation-produced tensor perturbations are stochastic in nature and char- 
acterized by a gaussian random variable. The Fourier components h]^ are 
drawn from a distribution whose statistical expectation is 

{Kh{) = PT(A:)(27r)«5(^)(k - q)5,,-; (19) 

where the gravity-wave power spectrum is defined as 



kr 



At = f (20) 
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The quantity V50 is the value of the scalar potential around 50 e-folds before 
the end of inflation when the modes of astrophysical interest, say A ~ 1 Mpc— 
lO^Mpc, were excited and crossed outside the horizon; see 0. [We note 
that the ensemble average and power spectrum are related by, = 
(27r)3Py(A;).] 

Using Eqs. (|18|, |1^) it is a simple matter to compute the energy density 
today in stochastic gravitational- wave background produced by inflation: 

where pcrit = SH^/SttG ^ 1.05/i^ x lO^eVcm"^ is the present value of the 



critical density. This agrees with previous results in the limit k <^ /ceq, 
and is shown in Fig. 3. 



The spectrum of gravitational waves is cut off at wavenumber |jT4| 

mp, llO"GeV; llO"GeV; * ' ^ ' 

where A^^ is the vacuum energy at the very end of inflation, Trh is the reheat 
temperature, and Tq = 2.726 K is the present temperature. The scale /cmax 
corresponds to the scale that crossed outside the horizon just as inflation 
ended, and is the highest frequency gravity wave produced. 

If 100% of the vacuum energy is converted into radiation at reheating, 
then inflation is followed immediately by the radiation-dominated epoch and 
Trh — M.- For "imperfect reheating," there is an epoch between the end of 
inflation and the beginning of radiation domination where the energy density 
of the Universe is dominated by the coherent oscillations of the scalar field 
responsible for inflation, and Trh < M.- In this case, the energy density in 
gravitational waves decreases as from k = k^, ^ Tq (Tmi/mpi) to k = fcmax 
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[0]; the scale crossed back inside the horizon just as the Universe became 
radiation dominated. 

The contribution of tensor perturbations to the variance of the multipoles, 
which arises solely due to the Sachs- Wolfe effect 0, ^, is given by 

(la^^P) = 36vrM^^^ T mu)\' \T{u/u^^)\' (23) 

1 £ — 1 Wo U 



where 



Fi{u) = - dy 



{u - y f 



(24) 



'('t.ss/to)" \ y 

u = kTo,y = kr, and Ueq = k^qTQ. 

This expression is identical to previous results with the exception of the 
inclusion of the transfer function to properly take account of short-wavelength 
modes, k > /ceq- Since \T{k/k-E(:i)\'^ ^ ^^r k > k^Q, without the trans- 
fer function the contribution from these modes has been underestimated. In 
Fig. 4 we show our results for the variance of the multipole moments, com- 
pared to those which do not take the short- wavelength modes into account 
properly. The correction to previous results is important even for values of / 
as small as 10 or so. 

The tensor contribution to the quadrupole CBR temperature anisotropy 
is given by 



^ oj^^a^ ^ 0.606^; (25) 

Jo-T 47r mpi^' ^ ^ 



V To J Q_T 47r mpi^ 

where the integrals in the previous expressions have been evaluated nu- 
merically; for the quadrupole term, including the transfer function does 
not make a significant difference. Using the COBE DMR measurement 
for the quadrupole anisotropy (derived for scale- invariant perturbations), 
ATq = 17 ± 4/iK as a rough upper limit to the tensor contribution. 
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Vso can be bounded |T^, 

V^so ^ 6.4 X lO^^^mpi^ ~ (3.5 x 10^*^ GeV)^ (26) 

This imphes that inflation takes place at a subPlanck energy scale (at least 
the last 50 or so e-folds that are relevant for us) and that Qqw can be at 
most about 10^^°. In models of first-order inflation another, a more potent 
source of gravity waves is produced by bubble collisions during reheating, 
though this radiation is peaked over a very narrow range of frequencies, 
k^2x 10^1 (Trh/10^^ GeV) Mpc"^ |17|. 



For reference, the Sachs- Wolfe result for the scalar modes, which is the 
dominant contribution to the scalar-mode produced anisotropy for I < 100, 
is 



ZTT JO U 



A 



10247r3 



50 



^'50 



AT\2 _ 5(|a2n^^ 327r V, 



To Jq-s 47r 45 mpi"^ x\ 



2 ' 
50 



(27) 



where u = krQ, ueq = /^eq, and Ts{y) is the transfer function for scalar 
perturbations, which depends upon the matter content of the Universe. For 
cold dark matter the transfer function is HTSl 



^ . ^ ln(l + 0.146y)/0.146y 

^^^^ [l + 0.242y + y2 + (o.340y)3+(0.417y)4]i/4- ^ ^ 

The scalar and tensor contributions to a given multipole are dominated 
by wavenumbers krQ ^ I. For scale-invariant perturbations and small /, both 
the scalar and tensor contributions to /^(|a;mP), are approximately constant. 
The contribution of scalar perturbations to /^(|a;mP) begins to decrease for 
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I = Ieq ~ 100 since the scalar contribution to these multipoles is dom- 
inated by modes that entered the horizon before matter domination and 
are suppressed by the (scalar) transfer function. The contribution of tensor 
modes to l'^{\aim\'^) begins to decrease for I ~ tq/tlss ~ Vl + -^lss (~ 35 for 
^hss ~ 1100) because the tensor contribution to these multipoles is domi- 
nated by modes that entered the horizon before last scattering (and hence 
decayed as until last scattering). The behaviour of the multipole am- 
plitudes is just as expected from the qualitative picture, cf. Fig. 1, and is 
illustrated in Fig. 4. 

2.3 Finite thickness of the last- scattering surface 

Last scattering is not an instantaneous event that occurs simultaneously ev- 
erywhere in the Universe. The last-scattering surface has a finite thickness, 
in comoving distance from the observer, ax, and in red shift, cr^. This fact 
leads to the damping of the contribution of modes with large wavenumber 
because the contribution to CBR anisotropy in a given direction averages 
over last-scattering events taking place over a finite distance, which washes 
out short-wavelength modes, k > ko — cr~^. This physics is not included 
in our analysis. Since the Ith multipole is dominated by the contribution of 
wavcnumbers /ctq ~ /, this damping only affects multipoles / > Tq/cTx, multi- 
poles that were already very small. Thus, our neglect of the finite thickness 
of the last-scattering surface does not affect our results in an important way. 

To be more specific, in the case of standard recombination, the differential 
probability that the last-scattering event occurred at comoving distance x — 
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the visibility function — can be approximated by a gaussian (see e.g., |TT 



^ ^ ^e-(-) ^ Cexp[-(x - a;Lss)V2^'], (29) 

where C ~ 182ifo; 3;lss — '^Hq^ is the center of the last-scattering surface, 
and ax = 2ifo^^/910 is the thickness of the last-scattering surface. (In red 
shift, the thickness is o"^ ~ 80.) Here s{x) is the optical depth from our 
position to a point comoving distance x from us, 

s{x) = j ne{x)aT — dx, (30) 

where ctt — 0.66 x 10~^^cm^ is the Thomson cross section and rig is the 
number density of free electrons (~ Xgn^). The damping scale associated 
with the thickness of the last-scattering surface is k^To = tq/cTx — 910, which 
leads to the suppression of the angular power spectrum for / > 910, where it 
is already very small, cf. Fig. 4. 

In the case of no, or only partial, recombination the last-scattering surface 
is very thick, ~ -^lss- The visibility function can be (less well) approxi- 
mated by a gaussian, 

dV 

— ~ Cexp[-(x - XLss)V2^^'], (31) 

where 

C = l/v^a,; XLss^2i7o"'[l-l/Vl + ^LSs]; a, ~ 2i/o" Vl2(l + ^lss)- 
The red shift of the last-scattering surface 

1 + zi^ss - {omx,nBh)'^/\ 

In the case of nonstandard recombination the damping scale kcTo ~ 4^zlss is 
smaller — because the last-scattering surface is thicker — but again, this damp- 
ing only affects multipoles that were very small anyway, / > 4:y/zi^ss ^ v^^lss, 
cf. Fig. 5. 
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2.4 Comparison with numerical results 

In Ref. [0] the coupled Boltzmann equations for the CBR intensity in a per- 
turbed FRW model were solved numerically to compute the CBR anisotropy 
that arises due to both scalar and tensor perturbations. The authors of 



Ref. have very gracious about comparing their preliminary results with 
us. Where comparison is possible, the agreement is always qualitatively very 
good, though there are quantitative disagreements: for values of / where the 
multipoles are of significant size at most about 20%. As we shall describe, 
we believe that we understand the reason for these disagreements. 

For very large /, corresponding to I > k^To, the numerical results of 
Ref. [|1^] fall off more rapidly than ours; as discussed above this is because 
we have not taken account of the damping due to the finite thickness of 
the last scattering surface. The discrepancy here — though large — is of little 
practical concern as the multipoles are very small for such values of /. 

A more important discrepancy arises in cases where last-scattering occurs 
before the Universe is "very matter dominated." The matter and radiation 
energy densities become equal at a red shift zeq — 2Ah'^ x 10^; for h = 0.5 
and standard recombination 2;lss and zeq differ only by a factor of 5.4. This 
means that the Universe is not well approximated as being matter dominated 
at last scattering. In particular, the growing mode gravity-wave perturbation 
is not given by 3ji{kT)/kT and 



TLBS = 2-^0 ^ (l + \/Req 



Rlss + Req — V -Req 



which differs from the matter-dominated result, tlss = '^Hq^ ^/ Ri^ss by a 
factor of 2/3. To correct Eq. ( P5D one would have to: (i) Use the correct 
expression for tlss; which is easy to do; (ii) Modify the transfer function. 
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taking into account that it is a function of red shift, which means that it can 
no longer be taken out of the inner integral; and (iii) Replace the conformal- 
time derivative of the growing-mode eigenfunction, the j2{y)/y term, by the 
proper expression which must be evaluated numerically and is a separate 
function of r and k (points (ii) and (iii) are not unrelated). Needless to say, 
these modifications eliminate the advantages of the analytic approach. 

The parameter that controls the size of the error made by assuming that 
the Universe is matter dominated at last-scattering is 2;lss/-2^eq- In the worst 
case considered, standard recombination and h = 0.5, this parameter is about 
0.2, which is about equal to the maximum difference between our analytical 
results and the numerical results of Ref. |T^. When zlss/-2eq is smaller, 
e.g., larger Hubble constant or nonstandard recombination, the discrepancy 
is much smaller. For example, with h = 0.5 and zlss = 76, the differences 



between the numerical results of Ref. |]T0| and our analytical results is only 
a few per cent (until the multipoles become very small due to the damp- 
ing associated with the finite thickness of the last-scattering surface). The 
advantage of the analytic approach is ease of calculation, particularly the 
ability to study a variety of scenarios. 



2.5 Generalized for tilt 

We now extend our results to the general case where the inflationary per- 
turbations are not precisely scale invariant. Since the ratio of the tensor to 
scalar contribution to the CBR quadrupole anisotropy increases with tilt, 
(kimP)/ (kfmP) ~ —7nT {riT = for scale invariance), the case of non scale- 
invariant perturbations is very relevant in discussing tensor perturbations. 
The deviations from scale invariance are most significant for steep potentials 
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(e.g., exponential potentials) and potentials whose steepness changes rapidly 
(e.g., low-order polynomial potentials). To lowest order in the deviation from 
scale invariance everything can be expressed in terms of the value of the po- 
tential V50, its steepness X50 = {rn-piV /V)5o, and the change in its steepness 



X' 



50' 



all evaluated about 50 e-folds before the end of inflation where the scales 



of astrophysical interest crossed outside the horizon during inflation fS], |T9|. 

Beginning with the tensor perturbations and using our previous notation, 
the power spectrum is given by 



Prik) 



A, 



riT 



kr 



8 V,o [r(f -|nT)]2 



3tt mpi^ (1 



|nr)2"T[r(|)]2"^^o ' 



X 



50 

Stt' 



(32) 



Here ^50 is the scale that crossed outside the horizon 50 e-folds before the 
end of inflation, rij- measures the deviation from scale invariance, and the 
expression for At includes the 0{nT) correction. (All formulas simplify if 
the potential and its derivatives are evaluated at the point where the present 
horizon scale crossed outside the horizon, i.e., ^50^0 ~ 1-) The variance of 
the multipole amplitudes is given by 



T{l + 2>) 

r(/-i) 



{\aim?) = ^Qti^ At 
where, as before, 

Fliu) = - I , dy 
u = kTo, y = kr, and ueq = kEQTQ. 



u 



riT 



\T{u/uE^)\^\Fi{u 



ji{u-y) 

{u - y f 



(33) 

u 



(34) 
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In Fig. 6 we show the scalar and tensor contributions to the angular 
power spectrum for spectra that are tilted by the same amount, n — l = nT = 
—0.15. This is an interesting case because the scalar and tensor contributions 
to the quadrupole anisotropy are essentially equal [|]. The effect of scale 
noninvariance is to tilt the angular power spectra, by approximately a factor 
of /""^ for scalar and /"^ for tensor. 

The analogous expressions for scalar perturbations are 

75H',+^ mpy^xl, [1 + Inr - i(n - l)]2-i[r(|)]2 ' 
^ - 1 ^50 , ^Pi4o . 

ilaimf) = A-^rl'- u-'\Ts{u/u^^)\'\3,{u)\'-- (35) 

ZTT JO U 

where u = krQ, ueq = k^q, and Ts{u/ueq) is the transfer function for scalar 
perturbations. 

By numerically integrating Eqs. (^,^) we can obtain expressions for 
the quadrupole anisotropy due to scalar and tensor perturbations in the non 
scale- invariant case: 



= 0.606-% (1 + 1.2n^) ; 
4:71 mpi'^ 



-7nT (1 + l.lriT + 0.1 [riT - {n - 1)]) . (36) 



(l«LP) 

We have taken /csoTq = 1, ^^lss = HOO (the results change very little for 
zlss = 76), and for the scalar case h = 0.5. Using these expressions, the fact 
that the scalar and tensor contributions to the quadrupole anisotropy add 
incoherently, and the COBE DMR quadrupole measurement, we can solve 
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for the variance of the tensor quadrupole, equivalently V^o/mpi^, in terms of 
the tensor tilt ut'- 



(l«LP) 



9.8 X 10-" 
1 -0.14n^^' 
1^50 6.4 X 10-" 



mpi^ l-0.14ny^' 



(37) 



These expressions indicate that the more tihed the gravity-wave spectrum 
is, the larger its amplitude is, as noted earlier 0]. 

The energy density in the stochastic background of inflation-produced 
gravitational waves today is given by 



dncw ^ , V50 [r(| - Inr)? 
dink mpi^ (1 - |nT)2"T[r(|)]2 

^ 1 _ il-i (1 + 0-^^^) ik/k,ornT{k/k^^)\' {kro)-'; (38) 

where the second expression follows from using the COBE DMR normal- 
ization to express the energy density in gravity waves in terms of the tilt 
parameter ut alone. In Fig. 3 we show the spectrum of stochastic gravi- 
tational waves for rix = —0.003, —0.03, —0.3, using the COBE quadrupole 
normalization. 

The total energy density in gravity waves increases with tilt (i.e., ut < 0), 
as does the tensor contribution to the CBR quadrupole anisotropy. However, 
this is not the entire story; the most sensitive "direct" probes of gravity 



waves, millisecond pulsars pO| and future Laser Interferometer Gravity-wave 
Observatories (LIGOs) , are only sensitive to gravity waves with very large 
wavenumber, kr^ = with ~ 26 for millisecond pulsars and ~ 48 for 
envisaged LIGO detectors. Because of the (fc/^so)"^ factor in Eq. ( P5| ) tilt 
depresses the energy density in gravity waves on the relevant scales. To be 
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more specific, for k ^ /ceq and h = 0.5 

dink nr-0.14 ^ ^ 

It is simple to sliow tliat tlie energy density in gravity waves on the scale 
k = c^Tq"^ is maximized for a value of rix ~ —1/N, at a value of about 
5 X 10^^^ /N. While the amount of tilt that maximizes the energy density in 
gravity waves on the scales relevant to both millisecond pulsars, ht — —0.04, 
and LIGOs, ut — —0.02, is quite reasonable in the context of inflationary 
models the predicted energy density in gravity waves is well below the 
sensitivity of either detector, about Qqw ~ 10~^° for advanced LIGO detec- 
tors and ~ 10^^ currently for millisecond pulsars. Thus, direct detection of 
the stochastic background of gravitational waves does not seem promising in 
the foreseeable future. 



3 Discussion 

The tensor contribution to the variance of the multipole amplitudes depends 
significantly upon the red shift of the last-scattering surface and less impor- 
tantly upon the red shift of matter-radiation equality. For scale-invariant 
gravity- wave perturbations /^(|a;mP) is roughly constant for / < V^T+^lss; 
then decreases as for / < y/1 + ^eq; and for / > a/T+^eq decreases as 
In the case of non scale-invariant tensor perturbations these results are 
modified by a factor of /"^. 

In Fig. 5 we show the tensor contribution to the angular power spectrum 
for no recombination and ^lss = 76. The dramatic fall off occurs at a rela- 
tively small value of /, around 10. Thus, the tensor angular power spectrum 
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can, in principle, be used to discriminate between standard recombination 
and no recombination, though due to cosmic variance, the finite "muhipole 
resolution" of experiments, and the difficulty of separating the scalar and 
tensor contributions to CBR anisotropy this is by no means a simple task. 
The angular power spectrum also depends upon the deviation of the ten- 
sor perturbations from scale invariance, both in its amplitude relative to the 
scalar perturbations and its dependence upon /. The angular power spectrum 
for tilted tensor perturbations is shown in Fig. 6. 

Finally, consider the cold dark matter + cosmological constant model 
(ACDM), proposed to reconcile a number of discrepancies of the CDM model 



with observational data p2[. It is characterized by: Q,b — 0.05, ficow — 0.15, 
Qq = ^Ib + ^coid = 0.2, Q/y ~ 0.8, and h ^ 0.8. For tensor perturba- 
tions most of the Sachs- Wolfe integral for the anisotropy arises near the last 
scattering surface where the effect of a cosmological constant is negligible. 
(This is not the case for scalar perturbations, and the formula for the Sachs- 
Wolfe contribution must be modiffed signiffcantly p3[.) Thus Eq. ( p3D for 
the contribution of tensor perturbations to the angular power spectrum is 
unchanged, with the following substitutions: 



Ho^Jno/REQ 



Req = 4.18 X 10"^ (l^o/i )"^ 3.27 X 10"^ 
To ^ H,' [' '^^ ^ 







tlss - 2^0 ' H^^/^l + zlss - 0.135i/o"\ (40) 

The angular power spectrum is shown in Fig. 7 for standard recombina- 
tion in the ACDM model. Unfortunately, the difference between the ACDM 
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and CDM models is not great. 

Finally, it is also straightforward to modify our results for the energy 
density in gravitational waves today for the ACDM model: 

(41) 

Since the main change is to reduce the energy density on a given scale by a 
factor of ~ 0.04, our conclusions about the direct detection of inflation- 
produced gravity waves remains. 
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FIGURE CAPTIONS 



Figure 1: The qualitative behaviour expected for CBR anisotropy arising 
from scalar and tensor perturbations. The horizon-crossing amplitudes of the 
scale and tensor perturbations are taken to be Sg oc A^-^"")/^ and St oc A~"'^/^ 
respectively (scale invariance corresponds to n — 1 = nr = 0) . 

Figure 2: The evolution of the transfer function for gravity- wave pertur- 
bations; from top to bottom, the transfer function computed at red shifts 
^ = 0, 30, 100, 300, 1000. 

Figure 3: The energy density of the stochastic background of inflation- 
produced gravitational waves, in the scale-invariant limit (broken curve, ar- 
bitrary normalization), and for ut — —0.003, —0.03, —0.3, normalized to the 
COBE DMR quadrupole (solid curves). 

Figure 4: The tensor contribution to the angular power spectrum nor- 
malized to the quadrupole, and for reference the scalar contribution nor- 
malized to the quadrupole (broken curve). Tensor results are shown with 

(lower curve) and without (upper curve) the transfer function; in all cases 
^Lss = 1100 and h = 0.5. (a) For / < 100; dotted curve shows tensor results 
with transfer function for h = 1.0. (b) For / < 1000. 

Figure 5: The tensor contribution to the angular power spectrum nor- 
malized to the quadrupole, for 2;lss = 76 (solid curve) and for comparison 
-^Lss = 1100 (broken curve); in both cases h = 0.5. 

Figure 6: Same as Fig. 4, except for non-scale invariant perturbations, 
n — 1 = Ut = —0.15. The results here correspond roughly to those in Fig. 4 
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^tilted" by a factor of {1/2)-' 



Figure 7: The tensor contribution to the angular power spectrum normal- 
ized to the quadrupole, for cold dark matter + cosmological constant with 
^^Lss = 1100, VLx = 0.8, and h = 0.8 (solid curve) and for comparison, 
^hss — 1100, flo, and h — 0.5 (broken curve). 
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